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SOLVABILITY AND LIMIT BICHARACTERISTICS 


NILS DENCKER 

1. Introduction 

We shall consider the solvability for a classical pseudodifferential operator P on a C°° 
manifold X. This means that P has an expansion p m + p m _ 1 + ... where pk G S^ om is 
homogeneous of degree k and p m = cr(P) is the principal symbol of the operator. The 
operator P is solvable at a compact set K C X if the equation 

(1.1) Pu = v 

has a local solution u G T>'(X) in a neighborhood of K for any v G C°°(X) in a set of 
finite codimension. We can also define the microlocal solvability at any compactly based 
cone K C T*X, see Definition 12.81 

A pseudodifferential operator is of principal type if the Hamilton vector field 

n 

(1-2) H p V (hjK).,. . - d Xj pd ii 

3 =1 

of the principal symbol p does not have the radial direction (£, d%) at p _1 (0), in partic¬ 
ular H p 7^ 0 then. By homogeneity H p is well defined on the cosphere bundle S*X = 
{ (x,£) G T*X : |£| = 1}, defined by some choice of Riemannean metric. For pseudodif¬ 
ferential operators of principal type, it is known m 0 that local solvability at a point is 
equivalent to condition ('I') which means that 

(1.3) lm(ap) does not change sign from — to + 

along the oriented bicharacteristics of Re(ap) 

for any 0 ^ a G C°°(T*M ). Oriented bicharacteristics are the positive flow-outs of the 
Hamilton vector held ffR e (a P ) 7^ 0 on Re(ap) = 0. Bicharacteristics of Re ap are also called 
semi-bicharacteristics of p. 

We shall consider the case when the principal symbol is real and vanishes of at least 
second order at an involutive manifold S 2 , thus P is not of principal type. For operators 
which are not of principal type, the values of the subprincipal symbol p m -i at S 2 becomes 
important. In the case when principal symbol p = £i£ 2 , Mendoza and Uhlman [5] proved 
that P was not solvable if the subprincipal symbol changed sign on the X\ or x 2 lines 
when = £2 = 0. Mendoza |6j generalized this to the case when the principal symbol 
vanishes on an involutive submanifold having an indefinite Hessian with rank equal to the 
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codimension of the manifold. The Hessian then gives well-defined limit bicharacteristis 
over the submanifold, and P is not solvable if the subprincipal symbol changes sign 
on any of these limit bicharacteristics. This corresponds to condition (P) on the limit 
bicharacteristics (no sign changes) because in this case one gets both directions when 
taking the limits. 

In this paper, we shall extend this result to more general pseudodifferential operators. 
As in the previous cases, the operator will have real principal symbol and we shall consider 
the limits of bicharacteristics at the set where the principal symbol vanishes of at least 
second order. The convergence shall be as smooth curves, then the limit bicharacteristic 
also is a smooth curve. We shall also need uniform bounds on the curvature of the 
characteristics at the bicharacteristics, but only along the tangent space of a Lagrangean 
submanifold of the characteristics, which we call a grazing Lagrangean space, see (12.51) . 
This gives uniform bounds on the linearization of the normalized Hamilton flow on the 
tangent space of this submanifold at the bicharacteristics. Our main result is Theorem l2.91 
which essentially says that under these conditions the operator is not solvable at the limit 
bicharacteristic if the quotient of the imaginary part of the subprincipal symbol with the 
norm of the Hamilton vector field switches sign from — to + on the bicharacteristics and 
becomes unbounded as they converge to the limit bicharacteristic. 

2. Statement of results 

Let the principal symbol p be real valued, £ = p - 1 (0) be the characteristics, and £ 2 
be the set of double characteristics, i.e., the points on £ where dp = 0. Let { Lj }^( 1 be 
a set of bicharacteristics of p on S*X so that Tj C £ \ £ 2 are uniformly bounded in C°° 
when parametrized on a uniformly bounded interval (for example with respect to the arc 
length). These bounds are defined with respect to some choice of Riemannean metric on 
T*X, but different choices of metric will only change the bounds with fixed constants. In 
particular, we have a uniform bound on the arc lengths: 

(2.1) \Tj\ < C Vj 

We also have that Tj = {7 j(t) : t G Ij} with = 1 and \I 3 \ < C, then | 7 ^(t)| < Ck 

for t G Ij and all j, k > 1. Let p = p/|Vp| then the normalized Hamilton vector field is 
equal to 

Hp = \H p \- ] H p onp-\0)\J: 2 
then Tj is uniformly bounded in C°° if 

(2.2) |P|Vp) < C k on Tj Vj, k 

where Vp is the gradient of p. Thus the normalized Hamilton vector field H p is uniformly 
bounded in C°° as a vector held over T. Observe that the bicharacteristics have a natural 
orientation given by the Hamilton vector held. Now the set of bicharacteristic curves 
{ Tj }“ =1 is uniformly bounded in C°° when parametrized with respect to the arc length, 
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and therefore it is a precompact set. Thus there exists a subsequence T ]k , k —> oo, that 
converge to a smooth curve T (possibly a point), called a limit bicharacteristic by the 
following definition. 

Definition 2.1. We say that a sequence of smooth curves Tj on a smooth manifold 
converges to a smooth limit curve T (possibly a point) if there exist parametrizations 
on uniformly bounded intervals that converge in C°°. If p G C°°(T*X) 1 then a smooth 
curve T C £ 2 f| ^*X is a limit bicharacteristic of p if there exists bicharacteristics Tj 
that converge to it. 

Naturally, this definition is invariant, and the set { T ? may have subsequences 
converging to several different limit bicharacteristics, which could be points. In fact, if 
Tj is parametrized with respect to the arc length on intervals I 3 such that \Ij\ —> 0 , then 
we find that Tj converges to a limit curve which is a point. Observe that if Tj converge 
to a limit bicharacteristic T, then ( 12 . ip and ( 12 . 2 p hold for Tj. 

Example 2.2. Let Tj be the curve parametrized by 

[ 0 , 1 ] 9 t 7 j(t) = (cos (jt), sin (jt))/ j 

Since | 7 '(f)| = 1 the curves are parametrized with respect to arc length, and we have 
that Tj —> (0,0) in C°, but not in C°° since I 7 "(t)| = j. If we parametrize T,- with 
x — jt £ [0, j] we hnd that T^ converge to (0,0) in C°° but not on uniformly bounded 
intervals. 

Example 2.3. Let P have real principal symbol p = w\ — a(w') in the coordinates 
(w!,w') e T* R n \0, where W\ = a(w') = 0 at £ 2 and k > 2. This case includes the cases 
where the operator is microhyperbolic, then k = 2 and a(w') vanishes of exactly second 
order at £ 2 - If k is even then in order to have limit bicharacteristics, it is necessary that 
a ^ 0 in a neighborhood of S 2 , and then £ is given locally by w\ = ±a(w') 1 ^ k . We hnd 
that (12.2p is satisfied if H J ~w\ = 0 and = 0(\ Va|) for any j when w\ = 0. 

But we shall also need a condition on the differential of the Hamilton vector held H p at 
the bicharacteristics along a Lagrangean space, which will give bounds on the curvature 
of the characteristics in these directions. I 11 the following, a section of Lagrangean spaces 
L over a bicharacteristic T will be a map 

T 3 w ^ L(w) C T W (T*X) 

such that L[w) is a Lagrangean subspace in T„,£, Via G T, where £ = p - 1 (0). If the 
section L is C l then it has tangent space TL C T l {Ty{T*X)). Since the Hamilton vector 
held H p (w ) is in T W T C T W (T*(X)), we hnd the linearization (or hrst order jet) of H p at 
w eT is in T Hp{w) (T w (T*X)). 
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Definition 2.4. For a bicharacteristic T of p we say that a C l section of Lagrangean 
spaces L over T is a section of grazing Lagrangean spaces of Y if L C TE and the 
linearization (or first order jet) of H p is in TL at T. 

Observe that since L(w) C T W T, is Lagrangean we find dp{w) | L( . , = 0 and H p {w ) G 
L(ic) when w G T. The linearization of H p {w ) is given by the second order Taylor 
expansion of p at w and since L(w) is Lagrangean we find that terms in that expansion 
that vanish on L{w) have Hamilton field parallel to L. Thus, the condition that the 
linearization of H p {w ) is in TL(w ) only depends on the restriction to L(w ) of the second 
order Taylor expansion of p at w. Since L(w) C T W T ,, we find that Definition 12.41 is 
invariant under multiplication of p by non-vanishing factors because p(w) = dp(w)\ L ^ = 
0. Thus we can replace H p {w) by the normalized Hamilton field H p in the definition. 

Example 2.5. Let L be given by T 9 w K> ( w,L(w )) with coordinates chosen so that 
L{w) = { (x,£) G R 2n : £ = A{w)x} where A[w) : R n (-)■ R n is linear. Then we can 
parametrize the tangent space of L(w o) with 

L'(w 0 ) : T W0 Y x L(w 0 ) 9 ( 5w , Sz ) (->• (tc 0 , L(w 0 ), Sw, ( 5x , A'(w 0 )(5w)5x) + 5z) 

and L is C 1 if A G C 1 (T,R n ). In these coordinates, the linearization of H p (wq) is given 
by 

H' p {w o) : T Hp{wo) (T wo (T*X )) 9 {Sx, Sz, SQ ^ H p (w 0 )(5z, 5() + d x ^H p (w 0 )(5x, SO 
where H p (w 0 )(5z, 5() = d i p{w 0 )5z - d x p(w 0 )5( G L(w 0 ). 

By Definition 12.41 we find that the linearization of H p gives an evolution equation for 
the section L, see Example 12.61 Choosing a Lagrangean subspace of T^ 0 E at wq G T then 
determines L along T, so L must be smooth. Actually, L is the tangent space at T of a 
Lagrangean submanifold of E, see (I4.22[) . 

Example 2.6. Let p = r—({A(t)x , x) + 2{B(t)x, 0 + (C(t)£, 0) /2 where (x, £) G T*R n , 
A(t), B(t ) and C(t) G C°° are real n x n matrices, where A(t) and C(t) are symmetric, 
and let T = { (£, 0, 0, ^ 0 ) : t G / }. Then H p = d t at T, 

P -1 (0) = {t= {A(t)x , x)/2 + {B(t)x, 0 + (C(t)O 0/2 } 

and the linearization of the Hamilton field H p at (t, 0, 0, ^o) is 

(2.3) d t + (A(t)y + B*(t)rj, d v ) - ( B(t)y + C(t) V , d y ) (;y , 77 ) G T*R n 

The linearization of the normalized Hamilton field H p is the same as (12.3ft . since \H p \ = 1 
modulo quadratic terms in ( y,r /). Since dp = dr at T, a C 1 section of Lagrangean spaces 
in T r E is for example given by 


£(f) = { ( S , y, 0 , £(£)?/) : (s, 2 /) G R n } 
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where E(t) G C 1 is symmetric with E(0) = 0, and by choosing linear symplectic coordi¬ 
nates ( y , p) we can obtain any such section on this form. By applying (12.3ft on p — E{t)y, 
which vanishes on L(t ), we obtain that L(t ) is a grazing Lagrangean space if 

(2.4) Eft) = A(t) + 2 Re B(t)E(t) + E(t)C(t)E(t) 

with Re F = (F + F*)/2, see (14.20ft . Then by uniqueness L(t) is constant in t if and only 
if Aft) = 0, observe that then Aft) = Hessp| L(4) . In general, Hessp| L(f) is given by the 
right hand side of (12.4ft . 


Observe that we may choose symplectic coordinates (t, x;r, £) so that r = p and the 
fiber of Lfw) is equal to { (s, y, 0, 0) : (s, y) G R n } at w G T = { ft, 0; 0, £o) : t G / }. But 
it is not clear that we can do that uniformly for a family of bicharacteristics { R, }, for that 
we need an additional condition. Now we assume that there exists a grazing Lagrangean 
space Lj of Tj, Vj, such that the normalized Hamilton vector field Hp satisfies 


(2.5) 


dHp(w ) | 


Lj(w) 


< C for w G Tj V j 


Since the mapping Tj 9 w ha Lj(w) 


is determined by the linearization of Hp on Lj, 


thus by dHp(w) | L , (w y condition (12.5ft implies that Tj 9 w ha Lj(w) is uniformly in C\ 


see Example 12.61 Observe that condition (12.2ft gives (12.5ft in the direction of T w Tj C 
Lj{w). Clearly condition (12. 5 p is invariant under changes of symplectic coordinates and 
multiplications with non-vanishing real factors. Now if 0 7^ u G C°° has values in R n and 
c 0 = u/\u\ is the normalization, then 


du = du/\u\ — (u, du)u/\u\ 3 


This gives that (12.5ft is equivalent to 


( 2 . 6 ) 


UdA/p\ L 


< C\Vp\ 


on Tj Vj 


where ni; — v — ( v , Vp)Vp/| Vp| 2 is the orthogonal projection along Vp. Condition (12.6[) 
gives a uniform bound on the curvature of level surface p -1 (0) in the directions given by 
Lj over Tj. Observe that the invariance of condition (12.6ft can be checked directly since 
d(ap ) = adp and d 2 (ap ) = ad 2 p + dadp = ad 2 p on Lj over Tj. 

The invariant subprincipal symbol p s will be important for the solvability of the op¬ 
erator. For the usual Kohn-Nirenberg quantization of pseudodifferential operators, the 
subprincipal symbol is equal to 


(2.7) 


Ps — Pm— 1 


2 i 




and for the Weyl quantization it is p m - 1 - 

Now for the principal symbol we shall denote 


0 < min \ H P \ — Kj —> 0 
Li 


( 2 . 8 ) 


] -> 00 
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and for the subprincipal symbol p s we shall assume the following condition 

(2.9) min / Imp s |.Hp| _1 ds/\ log Kj\ — > oo j —> oo 

J 

where the integration is along the natural orientation given by H p on Tj starting at some 

O 

point Wj £ Yj. This means that the integral must have a minimum less or equal to zero in 
the interior of Tj, and starting the integration at that minimum can only improve (12.9ft . 
Thus we may assume that the integral in (12. 9 p is non-negative. Observe that if (12.9(1 
holds then there must be a change of sign of Irn p s from — to + on Yj, and 

(2.10) max(—l) 1 * 11 Imp s /\H P \ | log/y,j —> oo j —> oo 

O 

for both signs. Observe that condition (12.9(1 is invariant under symplectic changes of 
coordinates and multiplication with elliptic pseudoclifferential operators, thus under con¬ 
jugation with elliptic Fourier integral operators. In fact, multiplication only changes the 
subprincipal symbol with the same non-vanishing factors as \H p \ and terms proportional 
to |Vp| = \H P \. Now by choosing symplectic coordinates (t,x;r,£) near a given point 
tco £ Tj so that p = ar near in 0 with a ^ 0, we obtain that d Xj d^p = 0 at Tj and 
d t d T p = d t a = <9 t |Vp| at Tj near w 0 . Thus, the second term in (12.7(1 only gives terms in 
condition (12.9(1 which are bounded by 

(2.11) J <9t|Vp|/|Vp| ds/ 1 log(«j)| < | log(|Vp|)|/| log(Kj-)| <1 j > 1 

since K-j < |Vp| on Yj. Here a < b (and b > a) means that a < Cb for some C > 0. Thus 
we obtain the following remark. 

Remark 2.7. We may replace the subprincipal symbol p s by p m _ 1 in (12.90 . since the 
difference is bounded as j —> oo. 

We shall study the microlocal solvability, which is given by the following definition. 
Recall that H l °Jj(X) is the set of distributions that are locally in the L 2 Sobolev space 

Definition 2.8. If K C S*X is a compact set, then we say that P is microlocally solvable 
at K if there exists an integer N so that for every / £ H l °X(X) there exists u £ U'(X) 
such that K H WF (Pu - f) = 0. 

Observe that solvability at a compact set M C X is equivalent to solvability at 
by PQ Theorem 26.4.2], and that solvability at a set implies solvability at a subset. Also, 
by Proposition 26.4.4 in [4] the microlocal solvability is invariant under conjugation by 
elliptic Fourier integral operators and multiplication by elliptic pseudodifferential opera¬ 
tors. The following is the main result of the paper. 

Theorem 2.9. Let P £ T™(A") have real principal symbol cr(P) = p, and subprincipal 
symbol p s . Let { Tj : be a family of bicharacteristic intervals of p in S*X so that (12.5p 
and ((2.9(1 hold. Then P is not microlocally solvable at any limit bicharacteristics of { Yj } .. 
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To prove Theorem 12.91 we shall use the following result. Let ||w||(fc) be the L 2 Sobolev 
norm of order k for u G C£° and P* the L 2 adjoint of P. 

Remark 2.10. If P is microlocally solvable at T C S*X, then Lemma 26.4.5 in [4] gives 
that for any Y d A" such that T C S*Y there exists an integer v and a pseudodifferential 
operator A so that WF(hL) D T = 0 and 

(2.12) IMI(-IV) < C(||P*-u||(^) + ||u||(_7V—n) + 11^11(0)) u G C£°(Y) 

where N is given by Definition 12.81 Since that definition also holds for larger N , we may 
take v = N > 0. 

We shall use Remark 12.101 in Section [7] to prove Theorem 12.91 by constructing approxi¬ 
mate local solutions to P*u = 0. We shall first prepare and get a microlocal normal form 
for the adjoint operator, which will be done in Section [4] Then we shall solve the eikonal 
equation in Section [5] and the transport equations in Section [61 

3. Examples 

Example 3.1. Let P have principal symbol p = Y\jPj which is a product of real sym¬ 
bols pj of principal type, such that pj = 0 on £ 2 , V j, and pj ^ pk on £ \ £ 2 when j ^ k. 
We find that if T C pjT^O) then 

\H p \ = \H Pk \ \pj\ = \H Pk \qk 
o+k 

where > 0 on T \ £ 2 close to £ 2 . Then p satisfies (12.211 and (12.51) for any section of 
Lagrangean spaces in Tp^ 1 (0) at T\£ 2 by the invariance, since Vp*, ^ 0 and d a pk = 0( 1), 
Va. A bicharacteristic T C £ \ £ 2 of p is a bicharacteristic for p*. for some k. Then if p s 
is the subprincipal symbol and (12.91) is satisfied for a sequence of bicharacteristics of pk 
converging in C°° to £ 2 , we obtain that the operator is not solvable at any limit of these 
bicharacteristics at £ 2 . 

Example 3.2. Assume that p(x,£) is real and vanishes of exactly order k > 2 at the 
involutive submanifold £ 2 = { £' = 0 }, £ = (£', £") G R m x R n_m , such that the localiza¬ 
tion 

V ^ °> f")*7 ° 

\a\=k 

is of principal type when p ^ 0. Then the bicharacteristics of p satisfies (12.21) and (12.51) 
with Lj = {£ = 0} at any point. In fact, \dpp(x, £)| = and d X: p'p(x,C.) = 0(\^'\ k ) 

so we obtain this since Hp = dppd x t + 0(|£'|) and cRVp = C , (|£ , | fc_1 ) when |£'| <C 1 and 
|£| = 1. The operator is not solvable if the imaginary part of the subprincipal symbol 
fmp s changes sign from — to + along a convergent sequence of bicharacteristics of the 
principal symbol and vanishes of at most order k — 2 at £ 2 . In particular we obtain the 
results of [5] and [6J. 
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Example 3.3. As in Example 12.31 let P have real principal symbol p = w\ — a{w') in 
the coordinates (wi,w r ) G T* R n \ 0, where W\ = a{w') = 0 at E 2 . We find that (12.51) is 
satisfied if dwA = 0 and dVal = C?(|Va|) when w\ = 0. 

Example 3.4. Let Q be a real and hyperbolic quadratic form on T*R". Then by j2j 
Theorem 1.4.6] we have the following normal form Qi(xfa)+Q 2 (y, rj) where the symplectic 
coordinates (x,y-fa,r)) G T*(R m x R r ), 

k m 

Qi(x, 0=^2 Viitf + ^ 2 ) + $ Vj > 0 v i 

i=l j=k +1 

is positive semidefinite and Q 2 (y,y) is either — fa, poViVi with p 0 ^ 0 or 2rjir} 2 ~ fa- To 
simplify the notation, we will assume /j,j = 1, \/ j. 

Now the Hamilton vector field HQ(x,yfa,rj) = Hgfaxfa) + Hg 2 (y,rj) so the flow of Hq 
is a direct sum of the flows of Hq 1 and Hq 2 . For Hq 1 it is a direct sum of the flows on 
the circles in [xjfaj) of radius rj —> 0, j < k, and the flows on x 3 lines, k < j < m. The 
circles can only converge in C°° to the origin if the radii r 3 go to zero fast enough, see 
Example 12.21 The possible limits are then points or lines in the Xk+i, ■ ■ ■, x m space. 

In the case Q 2 {y-, r i) — —yl we find that the limit bicharacteristics are given by 

(3.1) y\ = At, Xj = Xjt + dj, k < j < m 

where A 2 = J2'j‘=k+i We can only find a Lagrangean space satisfying (j2.5[) if Hj — 0, V j, 
since such a space cannot be contained in the subspace { Xj — Cj = 0 }. Then the grazing 
Lagrangean space can be chosen as { (x, y\ 0, 0) } at every point of the bicharacteristic. 
Theorem 12.91 then gives that the operator with homogeneous principal symbol equal to 
Q(x,y,£,v) no ^ solvable when the imaginary part of the subprincipal symbol changes 
sign on the lines on given by (13. lj) . which also follows from the results in |6J. 

If Q 2 {ypl) = yim, then <3 -1 (0) = {yiTji = -Qi(x,£) = -A 2 } where \H Ql \ = A. The 
Hamilton vector field of Q 2 is given by Hq 2 = y±d yi —rjid m which has flow 1 1 —> (i/ie*; 
with starting point (y\,r]i). Here yipi = —A 2 < 0 which gives |?/i| + |r/i| > lyiVi] 1 ^ 2 > 
We find that |Rq| = L — \yi\ + |t?i|, so in order for the bicharacteristic to converge in 
C°°, we find from (12.2ft that (\yi\ + |?/i|)/L fc = L l ~ k < 1, \/k. But then \yi\ + \r]i\ -fa 0 
so the bicharacteristics do not converge in C°° to a limit bicharacteristic at E 2 . Observe 
that a hyperbolic operator with Hessian of the principal symbol equal to Q(x,yfa,rj) 
at S 2 H S* R K is effectively hyperbolic and is solvable with any lower order terms, see the 
Notes to Chapter 23 in [4], 

Finally, when Q 2 (y, i/) = 2 rjir }2 ~ y\i then the characteristics in the (y,rj) variables are 
{yiy 2 = (y'i + A 2 )/2 } where A 2 = Q i. We find Hq 2 = 2p 2 d yi + 2r]id y2 + 2 y 2 d m so rji is 
constant on the orbits. Note that y\ + A 2 = 2r\\r\ 2 < fa + fa. Thus when faf > \rj 2 \ 
we find that |? 7 i| > \y 2 \ + A on £. Now |Rg 2 VQ 2 | > yfa so in order for (\2.2\i to hold 
we must have |?/i| > 1, so the characteristics will in that case not converge in C°° to 
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any limit bicharacteristic at T, 2 . When l?^) \vi\ we find that l^l 3> I 2 / 2 1 + A on £. A 
straightforward computation shows that applying Hq twice on VQ gives a factor 77 i/r] 2 . 
Thus, by (12.21) we find that the bicharacteristics converge in C°° to the y\ lines on £ 2 
only if |? 7 i| < 1 772 1 3 , which implies that |y 2 | + A < rfi. Then Hq = d yi + O(1 772 1) so the 
bicharacteristics can only converge to y\ lines at r \2 —> 0. An example is the case when 
y 2 — Vi — A = 0 and 772 —> 0, then Hq 2 = 2p 2 d yi and the grazing Lagrangean space can be 
chosen as { (x, y\ , 0; 0, 0, y 2 ) } at every point of the bicharacteristics. Thus Theorem 12.91 
gives that the operator with homogeneous principal symbol equal to Q{x,y,£,rj) when 
l£| 2 + M 2 = 1 is not solvable when the imaginary part of the subprincipal symbol changes 
sign on yi lines at £ 2 . 


4. The normal form 


First we shall put the adjoint operator P* on a normal form uniformly and microlocally 
near the bicharacteristics Tj converging to T. This will present some difficulties since we 
only have conditions at the bicharacteristics. By the invariance, we may multiply with 
an elliptic operator so that the order of P* is m = 1 and P* has the symbol expansion 
p + po + ..., where p is the principal symbol. As in Remark 12.71 we may assume that po 
is the subprincipal symbol. Observe that for the adjoint the signs in f)2.9[) are reversed 
and it changes to 

(4.1) max J Imp 0 |i7p| -1 ds/\ logKjj —>■ — 00 j —> 00 

where Kj given by (12.81) . By changing Wj to the maximum of the integral in (14.11) only 
improves the estimate so we may assume that 

(4.2) J lmp 0 /\H p \ ds < 0 on Tj 

with equality at Wj G Tj. Since V Imp 0 and VH P are bounded on S*X and \H p \ > Kj 
on Tj, we find that \H P \ and Imp 0 /|i7 p | only change with a fixed factor and a bounded 
term on an interval of length -C Kj on Tj. Therefore, we find that integrating lmp 0 /\H p \ 
over such intervals only gives bounded terms. Thus, we find from (12.101) that 


(4.3) \Tj\ » Kj 

and that condition ( 12 . 9 p holds on intervals of length = Kj at the endpoints of Tj. 

Now we choose 

(4.4) 1 < A, = K J 1/e » Kj = A 7 * 

for some e > 0 to be determined later. Then we may replace) log Kj\ with logAj in 
(12.9j) (12. lOj) . By choosing a subsequence and renumbering, we may assume by (12. 9 j) that 

(4.5) min / lmp 0 /\H p \ ds < — j log Xj 
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and that this holds on intervals of length = Kj at the endpoints of Tj. Next, we introduce 
the normalized principal and subprincipal symbols 

p = p/\H p \ and Ps = Po/\H p \ 

then we have that H p | r G C°° uniformly, \H p \ = 1 on Tj and dH p | L is uniformly bounded 
at Tj. We hnd that condition (14.5p becomes 

(4.6) min j Im p s ds < —jlogXj 

Observe that because of condition ( 12 . 101 ) we have that dTj has two components since 
Imp., has opposite sign there, so Tj is a uniformly embedded curve. 

In the following we shall consider a fixed bicharacteristic It, and suppress the index 
j, so that T = Tj, L = Lj and k = Kj = \~ £ . Observe that the preparation will be 
uniform in j. Now H p G C°° uniformly on T but not in a neighborhood. By (12.51) we 
may define the first jet of p at T uniformly. Since T G C°° uniformly, we can choose 
local coordinates uniformly so that T = { (i, 0)) : t G / C R} locally. In fact, take a 
local parametrization 7 (t) of T with respect to the arclcngth and choose the orthogonal 
space M to the tangent vector of T at a point wq with respect to some local Riemannean 
metric. Then R x M 3 (t,w) 1 —>• j(t) + w is uniformly bounded in C°° with uniformly 
bounded inverse near (t 0 , 0) giving local coordinates near T. We may then complete t to 
a symplectic coordinate system uniformly. We can define the first order Taylor term of p 
at T by 

(4.7) g(t, w) = d w p(t, 0) • w w = (x,t,£) 

which is uniformly bounded. This can be done locally, and by using a uniformly bounded 
partition of unity we obtain this in a fixed neighborhood of T. Going back to the original 
coordinates, we hnd that g G C°° uniformly near T and p — g = 0(d 2 ), but the error 
is not uniformly bounded. Here d is the homogeneous distance to T, i.e., the distance 
with respect to the homogeneous metric dt 2 + \dx\ 2 + ( dr 2 + |d£| 2 )/((r, £)) 2 . But by 
condition (12.51) we hnd that the second order derivatives of p along the Lagrangean 
space L at T are uniformly bounded. We shall use homogeneous coordinates, i.e., local 
coordinates which are normalized with respect to the homogeneous metric. 

By completing r — g in (14. 7 p to a uniformly bounded homogeneous symplectic coordi¬ 
nate system near T and conjugating with the corresponding uniformly bounded Fourier 
integral operator we may assume that 

(4.8) T = { (i, 0; 0, £ 0 ) : tel} 

for |£o| = 1, some hxed interval I 3 0, and p = r modulo second order terms at T. 
The second order terms are not uniformly bounded, but dVp \ L is uniformly bounded 
at T by (12. 5p . Since dp = dr on T we hnd that H p | = D t and we may obtain that 

L — { (t, x ; 0, 0) } at any given point at T by choosing linear symplectic coordinates (x, £). 
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Let 

(4.9) q(t, w) = |Vp(t, iu)| > k = X~ £ at T 

and extended so that q is homogeneous of degree 0, which is the norm of the homogeneous 
gradient of p. Recall that A is a parameter that depends on the bicharacteristic T. We 
shall change coordinates so that (t, w) = 0 corresponds to the point (0, 0; 0, £ 0 ) £ T given 
by (14. HD , and then localize in neighborhoods depending on A. We have |Vp| = 1 at L, 
higher derivatives are not uniformly bounded but can be handled by the using the metric 

g £ = ( dt 2 + |d«;| 2 )A 2£ 

and the symbol classes / G S(m,g £ ) defined by d a f = 0(mX^ £ ), Va. 

Proposition 4.1. We have that q is a weight for g £ , q G S(q,g £ ) andp(t,w) G S( A 1_£ ,p £ ) 
when |w| < cA _£ for some c > 0 when t G I. 

This gives p = qp G S(qX 1 ~ £ , g £ ) when |u>| < cA _£ . Observe that b G S}_ ££ if and only 
if b G S(X^,g £ ) in homogeneous coordinates when |£| > A > 1. In fact, in homogeneous 
coordinates 0 this means that dfb = (9(|£| /i+ l“l £ ) when \z\ = 1. 

Proof. We shall use the previously chosen coordinates (t, w ) so that T = { (t, 0) : t G / }. 
Now d 2 p = 0(1), q > X~ £ at T by (14.9p and 

(4.10) dq = Vp • ( dS7p)/q when q ^ 0 

which is uniformly bounded. We find that q(s,w) = q(t, 0) when |s — t\ + |iy| < cA _£ for 
small enough c > 0, so q is a weight for g £ there. This gives that \p(t,w)\ < q(t,w)X~ £ , 
\X7p(t,w)\ = q(t,w ) and \d a p\ < 1 < qX £ for |a| > 2, which gives p G S(qX~ £ ,g e ) when 
\w\ < cX~ £ . 

We find from (14. lOji that dq = 0(qX £ ) when |w| < cA _£ , since Vp G S(q,g £ ) in 
this domain. By induction over the order of differentiation we obtain from (14.101) that 
q G S(q,g £ ) when |w| < cA _£ , which gives the result. □ 

Next, we put Q(t,w ) = X £ p(tX~ £ ,wX~ £ ) when t G I £ — {tX £ : tel}. Then by 
Proposition 14.11 we find that Q G C°° uniformly when |w| < 1 and t G I £ , d r Q = 1 and 
\d t>x ^Q\ = 0 when w = 0 and t G I £ . Thus we find \d r Q\ ^ 0 for |w| < 1 and t G I £ . By 
using Taylor’s formula at T we can write Q(t , x] r, f) = r + h(t, x; r, ^) when |w| < 1 and 
t G I £ , where h — V/i| = 0 at w = 0. By using the Malgrange preparation theorem, we 
find 

t = a(t,w)(r + h(t,w)) + s(t,x,ff) |in| < 1 t G I £ 

where a and s G C°° uniformly, a ^ 0, and on T we have a = 1 and s = |Vs| =0. In fact, 
this can be done locally in t and by a uniform partition of unity for t G I £ . This gives 


(4.11) 


a(t, w)Q(t , w) = t — s(t, x,f,) 


< l tel E 
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In the original coordinates, we find that 

A e p(t, w) = a _1 (f A £ , wX £ )(tX £ — s(tX £ , xA £ , £A £ )) 

and thus 

(4.12) p(t,w) = b(t,w)(r-r(t,z,£)) |w| < A _e tel 

where 0 7 -be S(l,g £ ), r(t,x,£) = X~ £ s(tX £ , xX £ , £A £ ) e S(X~ £ ,g e ) when \w\ < A _£ , and 
t e /, b = 1 and r = |Vr| = 0 on T. By condition (12.51) we also find that 

(4.13) |dVr| L | < C when w — 0 and tel 

Recall that p = p/q , where q e S(q,g e ) when |iu| < A _£ . 

In the following, we shall denote by T the rays in T*X that goes through the bichar¬ 
acteristic. By homogeneity we obtain from (14.121) that 

b^q-'pit, x; t, f) = t - r{t, x, f) 

where 6” 1 e S®_ £ £ , q~~ x e Sf_ ££ and r — r e S\z e ££ when |£| > A and the homogeneous 
distance d to Y is less than c|£| _£ for some c > 0 . 

Next we take a cut-off function \ £ <S'i_ ee supported where d < |£| _£ and |£| ^ A so 
that b > Co > 0 in supp x and X = 1 when |£| > cA and d < c|£| _£ for any chosen c > 0. 
We let B = x& _ 1 g _1 e S{_ £ e and compose the pseudodifferential operator B with P*. 
Since P* e v h ) 0 we obtain an asymptotic expansion of BP* in S\tZ^ e 3( ' ] e> for j = 0, 1, 
2,... when d < |£| _£ . Bat actually the symbol is in a better class. The principal symbol 
is 

(r - r(t, x , 0)X e S{Zt, £ d < |£|" £ 

and the calculus gives that the subprincipal symbol is equal to 

(4.14) l -H v (xb~ l q~ l ) + Wo 

where p 0 is the subprincipal symbol of P*. As before, we shall use homogeneous co¬ 
ordinates when |£| > A. Then Proposition 14.11 gives p = qp e S(qX 1 ~ £ , g £ ) and since 
X' 6 _ 1 g _1 G S(q~ 1 ,g £ ) we find that (14.141) is in S(X £ ,g £ ) when d < A _£ and |£| > A. The 
value of H p at T is equal to qd t so the value of (14.141) is equal to 

(4-15) \-Aqlq + Po/q = when |f | > A at T 

2 1 2\X7 h p\ \X7 hP\ 

where |V*p| = y/\d x p\ 2 /|£| 2 + |%i | 2 is the homogeneous gradient, and the error of this 
approximation is bounded by A 2£ times the homogeneous distance d to T. In fact, since 
dp e S(q,g £ ) by Proposition 14.11 we have H p q~ l G S(X £ ,g e ) and po/q G S(q~ 1 ,g e ). 
Observe that Po/ 1V xp is equal to the normalized subprincipal symbol of P* on S*X. 
This preparation can only be done when |£| > A and the homogeneous distance to T is 
less than c|£| _£ , but we have to estimate the error terms. 
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Definition 4.2. For e < 1/2 and R G where g > e and 5 < 1 — e, we say that 
T*X 3 WF e(R) if the symbol of R is (9(|^| _Af ), V1V, when the homogeneous 

distance to the ray { ( x 0 , g£, 0 ) : g G R + } is less than c|£| _£ for some c > 0. 

By the calculus, this means that there exists A G ^i_ e £ so that A > c > 0 in a neighbor¬ 
hood of the ray such that AR G T~ 7V for any N. This neighborhood is in fact the points 
with fixed bounded homogeneous distance with respect to the metric g £ when A = |£|. 
It also follows from the calculus that this definition is invariant under composition with 
classical pseudodifferential operators and conjugation with elliptic homogeneous Fourier 
integral operators since the conjugated symbol is given by an asymptotic expansion by 
the conditions on g and 5. We also have that WF e (R) C WF(1?) when R G 

Now we can use the Malgrange division theorem in order to make the lower order terms 
independent on r when d < |£| _£ < A _£ , starting with the subprincipal symbol p 0 G S{_ £ E 
of BP* given by (14.141b Then rescaling as before so that Qo(t,w) = X~ £ po(tX~ £ , wX~ £ ) 
we obtain that 

Q 0 (t,w) = c(t,w)(r - s(t,x,£))+ q 0 (t,x,£) M<1 t G I £ 

where s is given by if 1.1 II) . and c and qo are uniformly in C°°. This can be done locally 
and by a partition of unity for t G I £ . We find in the original coordinates that 

(4.16) po(t,w) = c(t,w)(T-r(t,x,£)) + q 0 (t,x,£) d < X~ £ tel 

where q 0 {t,w) = X £ q 0 (tX £ ,wX £ ) G S(X £ ,g £ ) and c(t,w ) = X 2£ c(tX £ ,wX £ ) G S/A 2 '/^). By 
using a uniform partition of unity in T ( j ) _ £e , we obtain (14.16j) uniformly when |£| > A and 
the homogeneous distance to V is < |£| -£ with c G S 2 by homogeneity, q 0 G *S'[_ E£ and 
Qo — Po at T. Now the composition of the operators having symbols c and r — r gives 
error terms that are in I ££ when |^| > A. Thus if e < 1/3 then by multiplication with 
an pseudodifferential operator with symbol 1 — c we can make the subprincipal symbol 
independent of r. By iterating this procedure we can successively make any lower order 
terms independent of r when |£| > A and the homogeneous distance d to Y is less than 
c|£| _£ . By applying the cut-off function x we obtain the following result. 


Proposition 4.3. By conjugating with an elliptic homogeneous Fourier integral operator 
we may obtain that T is given by (14.81) . If 0 < e < 1/3 then for any c > 0 we may multiply 
with an homogeneous elliptic operator B G v b £ _ ee to obtain that BP* = Q + R G T}l ££ 
where rp|WF £ (l?) = 0 and the symbol of Q is equal to 


(4.17) r - r(t, x , 0 + q 0 (t, x, f) + r 0 (t, x, £) 

Here r G *S'/z £ £; q 0 G S £ l _ £ e and r 0 G S/ £ _~^ £ , r = 

p t \y h p(t,o)\ po(t, o) ( 2£ . 

2\Vhp(t, 0)| + \W h p(t,0)\ + [ > 


when d < c|^| e , |£| > cA and tel 
|Vr| = 0 on T, and q 0 is equal to 

when d < c|^| _£ , |^| > cA and tel 


(4.18) 
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where d is the homogeneous distance to T and \ Vhp\ = \f\d x p\ 2 /\£,\ 2 + |<%?| 2 is the homo¬ 
geneous gradient. This preparation is uniform with respect to A > 1. 

Observe that the integration of the term D t \VhP{t, 0) |/2| V^p(t, 0)| in (j4.18[) will give 
terms that are C?(log(| Vnp{t, 0)|)) = 0{\ log(A)| + 1), which do not affect condition fj2.9|) . 
We hnd from Proposition 14.31 that p(t,x; r,£) = r — r(t,x,f ) modulo terms vanishing of 
third order at T since r vanishes of second order at T. 

Recall that L is a smooth section of Lagrangean spaces L(w) C C T W (T* R n ), 
w E r, such that the linearization of the Hamilton vector held H p is in TL at T. By 
Proposition 14.31 we may assume that T = { ft, 0; 0, Co) : t E / }, 0 E /, and pft, x; r, C) = 
t — rft,x,f ) modulo terms vanishing of third order at T. Then we may parametrize 
L(t) = Lfw) where w = (i, 0, Co) for tel. Now since T* R n is a linear space, we may 
identify the fiber of T W (T* R n ) with T* R ra . Since Lfw) C T w Yj and w E T we hnd that 
r = 0 in Lfw). Since L(w) is Lagrangean, we hnd that t lines are in L(w). By choosing 
symplectic coordinates in (x,f) we obtain that L(0) = { (s, ?/;0,0) : (s,y) G R" }, then 
by condition (I2.5j) we hnd that <9;j:r(0, 0, £ 0 ) is uniformly bounded. Since r = 0 on L(t) 
and L(t) is Lagrangean we hnd for small t by continuity that 

(4.19) L(t) = { (s,y;0,A(t)y) : {s,y) e R n } 

where A(t) is real, continuous and symmetric for tel and A(0) = 0. Since the lineariza¬ 
tion of the Hamilton vector held H p at T is tangent to L, we hnd that L is parallel under 
the how of that linearization. Since L(t) is Lagrangean it is only the restriction of the 
second order jet of r(t,w ) to L(t) that determines the evolution of t i— L(t). For (I4.19P 
this restriction is given by the second order Taylor expansion of 

R(t, x) = r(t, x, + A(t)x) 

thus d 2 R(t,0) is uniformly bounded by condition (12.5^ . The linearized Hamilton vector 
held is 

d t + ( d 2 x R(t , 0)x, df) 

= d t + (( d%r{t , 0, Co) + 2 R e(d x d^r(t, 0, Co)4l) + Rcfr(t, 0, Co)4l) x, <9 5 ) 

where R e B = (B + B t )/2 is the symmetric part of B. Applying this on £ — A(t)x, which 
vanishes identically on L(t) for tel, we obtain that 

-A 1 (t) + (%r(t, 0, Co) + 2 R e(d x d^r(t, 0, Co)^(^)) + A(t)d%r(t, 0, Co )A(t) = 0 

which gives the evolution of L(t). The equation 

(4.20) Aft) = ^r(t,0,Co) + 2Re(a x %r(t,0,Co)A(t)) + A(t)a|r(t,0,Co)^W -4(0) = 0 

is locally uniquely solvable and the right-hand side is uniformly bounded as long as 
A is bounded. Observe that by uniqueness, Aft) = 0 if and only if cFr(t,0,Co) = 0, 
V t (see also Example 12.61) . But since (14.20P is non-linear, the solution could become 
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unbounded if dir ^ 0 and dir ^ 0 so that ||A(s)|| —* oo as s —» t± G I. This means 
that the angle between L(t) = { (s, y; 0, A(t)y) : (s,y) G R n } and the vertical space 
{ (s, 0; 0, r]) : (s, rj ) G R n } goes to zero, but this is just a coordinate singularity. 

In general, since we identify the fiber of T W (T* R n ) with T* R n we may define R(t, x, £) 
for each t so that 

(4.21) R(t,x,g) = r(t,x,& o + 0 when (0,x;0,£) G L(t) 

Then R = r on L and we find that 

(4.22) r- (R(t)z,z)/2 G C°° 

if z — (x, £) and R(t) = dlR(t, 0, 0) | L (i). Observe that we find from (12.51) that (14.221) 
is uniformly in C°° in z. We find that R(0) = d^r(t, 0,£o) and in general R(t) is 
given by the right hand side of (14.201) . Now we can complete t and (14.22ft to a uni¬ 
form homogeneous symplectic coordinates system so that T = { (i, 0,£ 0 ) : t G 1} and 
L(0) = { (s, ?/;0,0) : (s,y) G R n } . In fact, we may let x and £ have the same values 
when t — 0 and clearly H T is is not changed on T since then z — 0. This is a uniformly 
bounded linear symplectic transformation in (x, £) which is uniformly C 1 in t. It is 
given by a uniformly bounded elliptic Fourier integral operator F(t) which has uniformly 
bounded t derivative. We will call this type of Fourier integral operator a C 1 section of 
Fourier integral operators. This will give uniformly bounded terms when we conjugate 
a first order differential operator in t with F(t), for example the normal form which has 
symbol given by (14 .1 71) . For t close to 0 the section F(t) is given by multiplication with 
e (A(t)x,x)/ 2 , w } iere Afy solves (14.20ft . For general t we can put F(t) on this form after 
a linear symplectic transformation in (x, £). Observe that F(t) is continuous on local 
L 2 Sobolev spaces in x, uniformly in t, since it is continuous with respect to the norm 
||(1 + \x\ 2 + |.D| 2 ) fc w||, V k. In fact, it suffices to check this for the generators of the group 
of Fourier integral operators corresponding to linear symplectic transformations, which 
are given by the partial Fourier transforms, linear transformations in x and multiplication 
with e l ( Ax ’ x '> where A is real and symmetric. 

We find in the new coordinates that p — r — rr, where ri(t,x, £) is independent of r 
and satisfies dlri(t, 0, 0) \ L (t) = 0. This follows since 

p(t, T, £) = T - (R(t)z, z)/2 - n(t, X, r, £) 

where d T r\ = — {t,r±} = — {t, r } = 0 is invariant under changes of symplectic coor¬ 
dinates. Similarly we find that the lower order terms Pj(t,x,£) are independent of r 
for j < 0. Since the evolution of L is determined by the second order derivatives of the 
principal symbol along L by Example 12.61 we find that L(t) = { (t, x\ 0, 0) : (t, x) G R n }. 
Changing notation so that r = rr and p(t, x; r, £) = r — r(t, x, £) we obtain the following 
result. 
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Proposition 4.4. By conjugation with a uniformly bounded C 1 section of elliptic Fourier 
integral operators, corresponding to linear symplectic transformations in (x,£), we may 
assume that the symplectic coordinates are chosen so that the grazing Lagrangean space 
L{w ) = { (ft, x , 0, 0) : (ft, x ) G R n }, V w G T. which implies that d x r(ft, 0) = 0. 


We shall apply the adjoint P* of the operator on the form in Proposition 14.31 on ap¬ 
proximate solutions on the form 

M 

(4.23) u\(t, x) = exp(i X((x, £ 0 ) + w(t, x))) ^ <pj(t, x)X~ je 

3=0 


where |£ 0 | = 1 and the phase function u>(t,j G S(X~ 7e ,g 3e ) is real valued such that 
d x uj{t, 0) = 0 and cpj(t,x) G <5(1, < 7 $) has support where |x| < A” 5 . Here 5, e and g are 
positive constants to be determined later. The phase function u(ft,x ) will be constructed 
in Section O see Proposition 15.11 Observe that we have assumed that e < 1/3 in Propo¬ 
sition 14.31 but we shall impose further restrictions on £ later on. We shall assume that 
e + 5 < 1, then if p(t,x,£) G we obtain the formal expansion (see [7, Chapter VI, 

Theorem 3.1]) 


(4.24) p(t, x, D x )(exp(iX((x, f 0 ) + u(t, x)))<p(t, x)) 

~ exp(i X((x, f 0 ) + w(t, x))) ^ dfp(t, x, A (£ 0 + d x u(t, x)))lZ a {u, A, D)<p(t, x)/a\ 

a 

where U a (uj, X, D)<p(t, x) = D“(exp(iAw(t, x, y))<p(t, y))\ y=x with 

oj(t, x , y) = u(t, y ) - u(t, x) + (x - y)d x u(t , x) 

Observe that if \d x u\ -C 1 then this only involves the values of p(t,x,f ) where |£| > cA 
for some c > 0. Using this expansion we find that 


(4.25) P*(t, x , D)(exp(iX((x, £ 0 ) + u(t, x)))<p(t, x)) 

~ exp(iX((x,f 0 ) + u(t, x)))(Xd t ur(t, x) - r(t,x, A (£ 0 + d x u))<p(t,x) 

+ D t ip(t, x) - ^ d ^r(t, x, A (£ 0 + d x u))D x .(p(t , x) + q 0 (t, x, A (£ 0 + d x u))ip(t, x) 
j 

+ y X(f 0 + d x u))(D Xj D Xk ip(t,x) + iX(p(t,x)D Xj D Xk u(t,x))/2 + ...) 

jk 

which gives an expansion in S(X 1 ~ £ ~^ 1 ~ s ~ E \ gs), j > 0 , if 5 + e < 1 and e < 1/4. I 11 fact, 
since |£| = A every f derivative on terms in e gives a factor that is 0( A e_1 ) and every 
x derivative of p gives a factor that is 0(X s ). A factor A Dfu requires |a| number of f 
derivatives of a term in the expansion of P*, which gives a factor that is 0(A( 2 ~l a ^ 1_4£ ))_ 
Similarly, the expansion coming from terms in P* that have symbols in <Sf_ £e gives an 
expansion in SlZ J £ ^~ S ~ e \ j > 0. Thus, if 5 + e < 2/3 and e < 1/4 then the terms in the 
expansion have negative powers of A except the terms in (14.25I) . and for the last ones we 
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find that 

(4.26) ^ <%<%r(f, x -> Mfo + d x u))(D x .D Xk p + iXpD x .D Xk u) = £>(A 2<5+£_1 + X 3£ ~ s ) 

jk 

In fact, d^d^rft, x, A(£o + d x u)) = (9(A e_1 ) and D Xj D Xk uj = 0(X 2e d) when p ^ 0, since 
we have D Xj D Xk oj = 0 when x — 0, and d = 0( A _<5 ) in supp p. 

The error terms in 04.261) are of equal size if 25 + e — 1 = 3e — 5, i.e., 5 = (1 + 2e)/3. 
We then obtain 3 e — 5 = (7 £ — l)/3 < 0 if e < 1/7. Observe that in this case 1 — 5 — £ = 
(2 — 5e)/3 and 5 + e < 2/3 since e < 1/5. We also have that 1 — 4 e > 1 — 5 — £ since 
5 > 3 e. Thus we obtain the following result. 

Proposition 4.5. Assume thatuft,-) G S(X~ 7£ , g$ s ) is real valued and <9 2 a;(f,0) = 0, 
Pj(t,x ) G S(l,gs) has support where |x| < X~ s , for S, £ > 0. If 5 = (1 + 2e)/3 and 
£ < 1/7, then 04.25P has an expansion in S(X 1 ~ e ~^ 2 ~ 5e ^ 3 , g§), j > 0, and is equal to 

(4.27) exp(-iX((x, f 0 ) + w(t, x)))P*(t, x, D)(exp(iX({x, f 0 ) + u(t, x)))p(t, x)) 

~ (A d t ur(t, x) - r(t, x, X(f 0 + d x u)))p(t, x) 

+ D t p(t, x) - ^ d^rft, x, A(£ 0 + d x u))D Xj p(t, x) + q 0 (t, x, X(f 0 + d x cv))p(t, x) 

3 

modulu terms that are 0( A( 7£_1 )/ 3 ). 

In Section |6] we shall choose e — g — 1/10 which gives 5 = 2/5, 1 — 4 e = 3/5, 
1 - 5 - £ = 1/2 and (7e - l)/3 = -1/10. 


5. The eikonal equation 


The first term in the expansion 04.271) is the eikonal equation 


(5.1) 


d t co — s(t, x, f 0 + d x co) = 0 cu(0, x) — 0 


where s(t, x, £) 
tions: 

(5.2) 


A 1 r(t, x, A£). This we can solve by using the Hamilton-Jacobi equa¬ 
ls d t x = -<%s(t,x,£o + 0 
1 dtf = d x s(t,x,£ o + 0 


with (x(0),£(0)) = (x, 0), and letting d x u = f and d t uj = s(i, x, £ 0 + d x u>) withu;(0, x) = 0. 
Since s = Vs = 0 on T we find that d t x = d t £, = 0 when x = £ = 0, so by uniqueness 
d t u}(t, 0) = d x uj{t , 0) = 0. 

We shall solve the Hamilton-Jacobi equations by scaling. Recall that s(t, x, £) G 
S(X~ £ ,g e ) for some chosen 0 < £ < 1/3 in homogeneous coordinates by Proposition 14.31 
By Proposition 14.41 we may assume that L(t ) = { (i, x, 0, 0) }, V t, and d 2 s = 0 on T. Since 
s, ds and d 2 s vanish on T, Taylor’s formula gives 


(5.3) d$s(t,x,£ o + 0 = d x d^s(t,0,Co)x + <9fs(t, 0, £ 0 )£ + x, £)w, w) 
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where w = (x,f), d x d^s(t, 0,f o ) = 0(1) by (!23|) , <9|s(f, 0,£ 0 ) = 0(A £ ) and £1 G S(A 2e ,g £ ), 
since p(f, x;r, £) = r — r(t,x,£) modulo terms vanishing of third order at T. Similarly, 
we find 


(5.4) d x s(t, x, + 0 = d x d^s(t, 0, f 0 )f + (^(t, x, £)w, w) 


where g 2 £ S(A 2£ ,g £ ). 

Now we put (x, £) = 
transforms into 

(5.5) 


(yA 3e ,r/A 4e ). Then by using (15.3p and (15.4ft we find that (15.2(1 

( d t y = —B{t)y - C(t)g + afit, z) 

\ d t rj = B(t)rj + cr 2 (t,z) 


where z 
and 

(5.6) 


(y, rj), B{t ) = d x d^s(t, 0, £ 0 ) and C(t) = A £ d 2 s{t , 0, £ 0 ) are uniformly bounded, 

f cq/f, z) = A^ 3£ (g 1 (t,yA- 3£ 1 gA- 4£ )(y 1 A- £ g),(y,A- £ 7])} 

\a 2 (t,z) = -A- 2£ (g 2 (t,yA~ 3£ ,r]A- 4£ )(y, A~ £ r)), (y, A~ £ rj)) 


are uniformly bounded in C°° and vanish of second order in z. Then (15.5ft has a uniformly 
bounded C°° solution if 2 ( 0 ) is uniformly bounded. This means that if x(0) = 0( A _3e ) 
and £(0) = 0 then we find x = 0( A _3e ) and d x oj — £ = 0{ A _4e ) for any tel. The 
scaling gives that d x d x u> = 0(A0 4+3 I q D £ ) when |x| < A _3e . Since ca(t,0) e 0 we obtain 
that uj = 0(A _7e ) when |x| < A _3e , thus u>(t, •) G S(A~ 7e , g 3 £ ). 

Now by differentiating (15. ip twice we find that 


d t c%uj(t, 0) = 2 Re (d^d x s(t, 0, £ Q )d$.u(t, 0)) + c%u(t, 0)<9|s(f, 0, £ 0 )d%.u(t, 0) 

because d x u>(t, 0) = d^s(t,0, £0) = d 2 s(t, 0, £0) = 0. Since d 2 u(0,x) e 0 we find by 
uniqueness that d 2 u(t, 0) = 0. This gives that d t c 0 = 0(A _7e ) when |x| < A _3e , and we 
obtain the following result. 


Proposition 5.1. Let 0 < e < 1/3, and assume that Propositions |^.3| and |^./| hold. 
Then there exists a real G S(A~ 7s , g^fi satisfying d t u = A _1 r(f, x, A(£o + d x u )) 

when |x| < A~ 3e and tel such that co(t,0) = 0 and d 2 cu(t, 0) = 0. We find that the 
values of (t, x; Ad t co(t, x), A(£ 0 +d x cu(t, x))) have homogeneous distance to the rays through 
T which is < A _3e when Ixl < A _3e and tel. 


6. The transport equations 

The next term in (|4.27(1 is the transport equation, which is equal to 

(6.1) D p ip 0 + q 0 <p 0 = 0 at T 

where D p = D t - JN d^r(t, x, A(£ 0 + d x u(t, x)))D Xj = D t when x = 0 and 

(6.2) q 0 (t) = DfiVplfi, 0, £o)|/2|Vp(t, 0, £ 0 )| + PoM, £ 0 )/|Vp(t, 0, £ 0 )| = 0(A £ ) 
modulo 0( A 2e |x|) when |x| < A” 5 by (14.181) . Here uo(t,x) is given by Proposition 15. 11 
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Lemma 6.1. We have that 

D p = D t + y^ / (a ] (t),x)D Xi + R(t, x, D x ) 

3 

where R n_1 9 aj(t) = 0( 1) and R(t, x, D x ) is a first order differential operator in x with 
coefficients that are 0( A 3e |x| 2 ) when \x\ < A _3e . 

Proof. As before, we shall put s(t,x,£) = A _1 r(f, x, A£) G S(X~ £ ,g e ). Since <9 2 cu(f,0) = 0 
by Proposition 15.11 we have from Taylor’s formula that aj(t) = —d x d^s(t, 0, £ 0 ) which is 
uniformly bounded by (14.131) and Proposition 14.41 The coefficients of the error term R 
are given by the second order x derivatives of the coefficients of D p which are 

dj/fis + 2 Re {d x 8lsdloj) + d 2 x ud 3 sd 2 x oj + dfs<9 3 u; = 0(\ 3e ) 

when |x| < A _3e -C X~ £ by Propositions 14.31 and 15.11 which proves the result. □ 

We obtain new variables y in R"^ 1 by solving 

dtVj = Wj(t),y) yfi 0) = Xj Vj 

Then D t + (aj (t ), x)D Xj is changed into D t but D Xj = D yj is unchanged, and we will for 
simplicity keep the notation (t,x). The change of variables is uniformly bounded since 
a, = 0(1), so it preserves the neighborhoods |x| < X~ u and symbol classes S(X fl ,g u ), 
V fi, u. We shall then solve the approximate transport equation 

(6.3) D t ip 0 + q 0 (t)(po = 0 

where <^ 0 (0, a;) G S(l,gs) is supported where |x| < A -5 and q 0 (t) is given by (j6.2j) . If 
A -5 <C A _3e then by Lemma [6.1 1 the approximation errors will be in S(X 3£ ~ s , gf), so we 
will assume S > 3e. In fact, since d x maps S(l,g$) into S(X s ,g$ ) and \x\ < X~ 6 , we 
find R(t,x, D x )(po G S(X 3e ~ s , gf). If we put 5 = 4e then the approximation errors in the 
transport equation will be 0(X~ £ ). 

If we choose the initial data </? o (0,x) = <j>o( x ) — ^(A* 5 #), where <p G Cq° satisfies 
<^(0) = 1, we obtain the solution 

(6.4) ip 0 (t, x) = 4> 0 (x) exp (~iB(t)) 

where B' = q 0 and B(0) = 0. By condition (14. 2 p we find that \ip 0 (t,x)\ < |(^(A‘ 5 x)|, 
and |x| < A -8 in supp</?o, which also holds in the original x coordinates. Observe that 
D x ip o = (P(A' 5 I“I), Va, and we have from the transport equation that D t ip 0 = — q 0 (p 0 = 
0(X £ ) by (j6.2jl . Since D k q 0 = 0(X e( ~ k+1 ^) by Proposition 14.31 we find by induction that 
<Po G 5(1, gs). 

After solving the eikonal equation and the approximate transport equation, we find 
from Proposition 14.51 that the terms in the expansion (14.271) are 0(X^ £ ~ 1 ^ 3 ) + 0(X~ £ ), if 
e < 1/7 and 5 = (1 + 2er)/3 = 4e, and all the terms contain the factor exp(—iB(t)). We 
take e = 1/10, 6 = 2/5 which gives (7e — l)/3 = —e = —1/10. Then the expansion in 
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Proposition 14.51 is in multiples of A -1 / 2 , but since the terms of (14.27p are 0( A” 1 / 10 ) we 
will take g = 1/10. 

Thus the approximate transport equation for p\ is 

(6.5) D t tpi + qo(t)p i = A 1//10 i?i exp(—iB(t)) at F 

where R\ is uniformly bounded in the symbol class S/A -1 / 10 , ( 72 / 5 ) and supported where 
|x| < A -2 / 5 . In fact, R\ contains both the error terms from the transport equation (16.11) 
for p 0 and the terms that are 0(A -1 / 10 ) in (14.271) . By putting 

pi(t,x) — exp(— iB{t))(j)i(t, x) 

the transport equation reduces to solving 

(6.6) D t <j> 1 = X 1/W R 1 

with initial values 0i(O, a;) = 0. Then 0! G 5/1, .c/ 2 / 5 ) will have support where |x| < A -2 / 5 . 

Similarly, the general term in the expansion is pkX~ k ^ 10 where pk will solve the ap¬ 
proximate transport equation 

(6.7) D t p k + qo(t)<pk = \ k/w Rk exp(iB(t)) k> 1 

with R k is uniformly bounded in the symbol class 5/A _fc / 10 , ( 72 / 5 ) and is supported where 
|x| < A^ 2 / 5 . I 11 fact, Rk contains the error terms from the transport equation (16.11) 
for ipk-i and also the terms that are 0{ A _fc//1 °) in (14.27ft . Taking p k = exp {—iB{t))4>k we 
obtain the equation 

(6.8) D t (j) k = X k/w R k G S(l,g 2 / 5 ) 

with initial values 0*00, x) = 0, which can be solved with 0*, G 5/1, . 92 / 5 ) uniformly 
having support where |x| < A -2 / 5 . Proceeding we obtain an solution modulo 0(X~ N ^ 10 ) 
for any N. 

Proposition 6.2. Choosing 6 = 2/5, e = 1/10 and g = 1/10 we can solve the transport 
equations (16.1ft and (16.71) with p k £ 5/1, $ 2 / 5 ) uniformly having support where |x| < A -2 / 5 
and |i| < 1, W k, such that po(0, 0) = 1 and pk(0, x) = 0, k > 1. 

Now, we get localization in x from the initial values and the transport equation. To get 
localization in t we use that ImB/) < 0. Then we find that Re(— iB) < 0 with equality 
at t — 0. Near <9r we may assume that R e(—iB(t)) -C — log A in an interval of length 
0(X~ £ ) = 0( A -1 / 10 ) by (14.31) . (14.41) and (14.6ft . Thus by applying a cut-off function 

X (t)eS(l,A 1 / 5 dt 2 )cS(l,g 2/5 ) 

such that y(0) = 1 and y/f) is supported where (14.51) holds, i.e., where p k = 0(X~ N ), 
V k, we obtain a solution modulo 0(X~ N ) for any N. In fact, if u\ is defined by (14.231) 
and Q by Proposition 14.31 then 


Qxu a = xQu a + [Q, x]u\ 
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where [Q, x] = DtX is supported where u\ = 0(X~ N ) which gives terms that are (D(X~ N ), 
V N. Thus, by solving the eikonal equation (15.ip for u> and the transport equations (16.71) 
for p k for k < 10 IV, we obtain that Qx u a = 0{X~ N ) for any N and we get the following 
remark. 

Remark 6.3. In Proposition 16.21 we may assume that p k (t,x) = 4> k (X 1 ^ 10 t, A 2 / 5 x) G 
S(l,g 2 /f), k > 0, where <f k G Cff has support where |x| < 1 and |i| < A 1//10 < A 2//5 , 
A > 1. 


7. The proof of Theorem 12.91 

For the proof we will need the following modification of [4] Lemma 26.4.14], Recall 
that V' r = {u G V : WF(u) C T } for T C T*R n , and that ||n||(fc) is the L 2 Sobolev norm 
of order k of u G . 

Lemma 7.1. Let 

M 

(7.1) u x {x) = A (n_1)<5/2 exp(iAMA £ a;)) <pj(X 5 x)X~ jK 

j =o 

with u G C°°(R") satisfying Imw > 0 and \du\ > c > 0, G C'“(R n ), A > 1 , e, 8, k 
and g are positive such that e < 5 < e + g. Here to and <pj may depend on A but uniformly, 
and tpj has fixed compact support in all but one of the variables, for which the support is 
bounded by CX s . Then for any integer N we have 

(7.2) ||ma|| ( -jv) < C\~ N ^ e+e ^ 

If <p o(x 0 ) 0 an d Imu(x 0 ) = 0 for some x 0 then there exists c > 0 and A 0 > 1 so that 

(7.3) IKH(-iv) > cX- {N+ ^ £+Q)+{n ~ 1)5/2 A > A 0 
Let S = Ha>i U, supp <Pj(\ •) and let T be the cone generated by 

(7.4) { (x,dw(x)), x G E, Imw(s) =0} 

then for any real m we find \ m u x -A 0 in V r so X m Au\ —> 0 in C°° if A is a pseudo- 
differential operator such that WF(A) fl T = 0. The estimates are uniform if oj G C°° 
uniformly with fixed lower bound on |dReo;|, and >p 3 G C°° uniformly. 

We shall use Lemma [7. II for u x in (14.231) . then oj will be real valued and T in (17.41) will 
be the bicharacteristic Tj converging to a limit bicharacteristic. 

Proof of Lemma \ 7. 1\ We shall adapt the proof of J4J Lemma 26.4.14] to this case. By 
making the change of variables y = X e x we find that 

M r 

MO = X in ~ 1)5/2 ~ n£ J2 x ~ jK / e^Avi-M^MX s ~ e y)dy 

3=0 •' 


(7.5) 
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Let U be a neighborhood of the projection on the second component of the set in (17.41) . 
When £/A £+e ^ U then for A 1 we have that 

U supp <Pj{ A 5 ' 6 -) 3y^ (A e 0 j(y) - (y, ^/A £ ))/(A £) + |£|/A £ ) 

3 

= Mi,)-{ 9 ,5/V + «))/(1 + |5|/A' + «) 

is in a compact set of functions with non-negative imaginary part with a fixed lower 
bound on the gradient of the real part. Thus, by integrating by part in (17.51) we find for 
any positive integer m that 

(7.6) |h A (0| < C m \-("-W2+ m (s-e)( X e + |^|/A £ )- m £/A £+e <£U A > 1 

This gives any negative power of A for m large enough since 5 < £ + g. If V is bounded 
and 0 ^ V then since U\ is uniformly bounded in L 2 we find 

[ Mora+i t\ 2 r N <% < cvt-™ r> i 

JtV 

Using this estimate with r = A £+e together with the estimate (17.61) we obtain (17.21) . If 
X € Cq° then we may apply (17.6[) to x u \, thus we find for any positive integer j that 

|x«a(OI < QA- (n - 1)5/2+i(5 - £) (A 0 + |£|/A £ )^ (GIL A » 1 

if W is any closed cone with T p|(supp x x IU) = 0. Thus we find that A m u\ —> 0 in V' T 
for every m. To prove (17.3[) we assume xq = 0 and take -0 G Cq°. If lmu;(0) = 0 and 
99 ( 0 ) ^ 0 we hnd 

X n(e+e)-(n-l)6/2 e -i\e Keu(O) ^ ^\ e+e -)) 

= f e iXe ^ x/xe) -^ 0)) ^{x)J2^j( x / xe+e ~ S ) X ~ jKdx 

j 

-G J e dKea ^ x) il3{x)^{{))dx A ^+00 

which is not equal to zero for some suitable 0 G C£°. In fact, we have tpj(x/\ £+e ~ s ) = 
<Pj( 0) + 0(\ 5 ~ £ ~ q ) —y <pj( 0) when A —> 00 , because 6 < e + g. Since 

II^(a £+ "-)II(tv) < c\ {N ~ n/2) ( £+e) 

we obtain that 0 < c < A( JV+ t)( £ +^) _ ( n_1 ) <5 / 2 ||' lt ||(_ JV ) which gives (17. 3 p and the lemma. □ 

Proof of Theorem\ 1 2. .91 Assume that T is a limit bicharacteristic of P. We are going to 
show that (12.121) does not hold for any is, N and any pseudodifferential operator A such 
that T D WF(A) = 0. This means that there exists approximate solutions 0 7 - Uj £ C^° 
to P*Uj = 0 such that 

(7-7) ||^||(-iv)/(||P%'|Im + IMI(-iv-n) + ll^ll(o)) 00 when j 00 


which will contradict the local solvability of P at T by Remark 12.101 
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Let Tj be a sequence of bicharacteristics of p that converges to T C S 2 and A j given 
by (14.3p and (j4.4[) with e to be determined later. Now the conditions and conclusions 
are invariant under symplectic changes of homogeneous coordinates and multiplication 
by elliptic pseudodifferential operators. Thus by Proposition 14.31 we may assume that the 
coordinates are chosen so that Tj = I x (0, 0, £j) with |£j| = 1, and for any 0 < e < 1/3 
and c > 0 we can write P* = Q + R where Tj D WF e (R) = 0 and Q has symbol 

(7.8) r - r(t, x, g) + q 0 (t, x, f) + r 0 (t, x, g) 

when |£| > c\j and the homogeneous distance to Tj is less than c|£| _£ . We have that 
R 6 Spp, t o G Sf I e 1 e , q 0 G Sf_ £ £ is given by (14.181) . and r G <5q_ £ e vanishes of second 
order at T 3 . 

Now, we may replace the norms ||m||( s ) in (17.71) by the norms 

IMI* = \\(d x ) s u\\ 2 = I (0 2a K r »0l 2 drd £ 

In fact, the quotient (£)/((t, £)) — 1 when |r| < |^|, thus in a conical neighborhood of T. 
So replacing the norms in the estimate (17.71) only changes the constant and the operator 
A in the estimate (12.121) . By using Proposition 14.41 we may assume that the grazing 
Lagrangean space Lj(w ) = { (s,y;0,0) : (s,y) G R n }, Vw G Tj, after conjugation with 
a uniformly bounded C 1 section F(t) of Fourier integral operators, then d^r = 0 at Tj 
and T D WF e (4) = 0. Observe that for each t we find that F(t) is uniformly continuous 
in local H s spaces, which we may use in (17.7|) after changing A. Also the conjugation 
of F(t) with the operator with symbol (I7.8f) gives a uniformly bounded expansion. By 
changing A again, we may then replace the local ||w|| s norms by the norms ||w||( s ) in (17.71) 
so that we can use Lemma 17.11 

Now, by choosing S = 2/5, £ = 1/10 and p = 1/10 and using Propositions 14.51 
15.1116.21 and Remark 16.31 we can for each r ; construct approximate solution U\ :J on the 
form (14.23ft so that Qu\ j = 0(X k ), for any k. The real valued phase function is equal 
to is (x,g) +ujj(t,x ) where \g\ = 1 and uJj{t,x) G <S'('/ 10 , g 3 / 10 ) and the values of 
(t, x; Xd t u>j(t, x), X(g + d x Uj(t,x))) have homogeneous distance to the rays through Tj 
which is < A~ 2 / 5 when |x| < Aj~ 2 ^ 5 , i.e., in supp u\ ) . Observe that if A 1 then we have 
that |f 0 + d x cu(t,x)\ > c in supip u\^ for some c > 0. We find that 

Idj (t, x) = Xj~ 7 ^ 10 iJjXj 3 ^ W t, Xj 3 ^ 10 x) 

where uj 3 G C°° uniformly so d x u = 0(XJ 2 ^ 5 ) when x = 0(XJ 2 ^ 5 ) and 

+Vj(t,x)) = X j 7/10 (X j 3/10 x, g) + X~ 4/10 uj j (X j 3/w t, Xj 3/w x) 

Thus 6 = 2/5, £ = 3/10, a = 1/10 and g = 7/10 in (17. ip so we find £ + £> = l><5>£. 

The amplitude functions ipkj{t,x) = /'fcj(A J 2 ' / % A ? 2 ^ 5 t) where 0^- G C^° uniformly 
in j with fixed compact support in x, but in t the support is bounded by CXj 2 ^ 5 , so u\. 
will satisfy the conditions in Lemma ITTl uniformly. Clearly differentiation of Qu\ j can at 
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most give a factor A j since e + q = 1 and S < 1. Because of the bound on the support of 
u\. we may obtain that 

(7-9) WQuxjWw = 0(\- N ~ n ) 

for the chosen u. 

If WF(A)fir = 0, then we find WF(A)f'|r j - = 0 for large j, so Lemma [7TI gives 
II^Ajll(o) = 0(\J N ~ n ) when j —> oo. Since x = 0(X J 2/5 ) in supp u \ :j , we hnd that 
the values of (t, x; Xd t tUj(t,x), X(£j + d x Uj(t,x))) have homogeneous distance to the rays 
through Tj which is < A -2 / 5 for x G supp u \., and this converges to T. Thus, if R G 
‘Sg/ioVio suc ^ that WF 1/10 (i?) (jTj = 0 then we hnd from the expansion (14.241) that all 
the terms of Ru\ vanish for large enough A,-. In fact, since X~ 2 ^ -C A7 1 ^ 1 ° for j 1, we 
hnd for any a and K that 

d a R(t,x-, Xj((0,Q + dt, x u}j(t,x))) = 0(Xj K ) 

in Uk SU PP Vk,j- As before, we hnd that \\Ru\ j ||(^) = 0(Xj N ~ n ) by the bound on the 
support of n a, so we obtain from (17.9ft that 

(7.io) \\P* u Xj\\n = o(xj N - n ) 

for the chosen u. 

Since £ + g = 1 we also hnd from Lemma 17.11 that 

\—N \—N(£+q) \ I. || s, . —(N+S.)(e+e)+(n—1)5/2 ,-V— n/2+(n— 1)/5 s. \—N—n /2 

Xj -Xj > || wa 3 .||(- jv ) ^ Xj — Xj — Xj 

when Xj > 1. We obtain that (17.71) holds for Uj = u\ j when j —> oo, so Remark 12.101 
gives that P is not solvable at the limit bicharacteristic T. □ 
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